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Abstract 

Rather than regarding the restriction of current lattice QCD simulations to 
quark masses that are 5-10 times larger than those observed, we note that 
fC) \ this presents a wonderful opportunity to deepen our understanding of QCD. 

Just as it has been possible to learn a great deal about QCD by treating N c 
as a variable, so the study of hadron properties as a function of quark mass is 
leading us to a much deeper appreciation of hadron structure. As examples we 
\ cite recent progress in using the chiral properties of QCD to connect hadron 
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masses, magnetic moments, charge radii and structure functions calculated at 
large quark masses within lattice QCD with the values observed physically. 



O 

£ . I. INTRODUCTION 



In striving to understand the properties of QCD the generalization to an arbitrary num- 



ber of colours, N c , particularly the limit N c — > oo (or "large iV c ") has been extremely 
valuable. It has even proven possible to distinguish between models of hadron structure and 
to guide the further developments of such models on the basis of their large N c behaviour 
0. Until recently it has generally been regarded as an unfortunate liability that current 
limitations on computer power restrict lattice QCD simulations with dynamical fermions to 
large quark masses. We would like to present a rather different view concerning the lattice 
data at large quark masses. In particular, we argue that like the behaviour as a function of 
N c , lattice results as a function of quark mass offer extremely valuable new insights into the 
nature of QCD and especially into hadron structure. 

To be a little more quantitative, the restriction to large quark masses in lattice simu- 
lations means typically 50 MeV or higher. Thus, in order to compare hadron properties 
calculated on the lattice one has to extrapolate as a function of quark mass (on top of all 
the other extrapolations, lattice spacing, lattice size, etc.) all the way to the physical light 
quark masses, around 5 or 6 MeV. Such extrapolations are complicated enormously by the 
fact that chiral symmetry is spontaneously broken in QCD. The mass of the pion, which is 
the Goldstone boson corresponding to this broken symmetry |§, behaves as: 
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(with m = m u = rrid 7^ 0), 



(1) 



as the quark mass, fh, moves away from zero - this is the Gell Mann-Oakes-Renner (GOR) 
relation. While Eq.(l) is, in principle, only guaranteed for quark masses, near zero, explicit 
lattice calculations show that it holds over an enormous range, as high as ~ lGeV. 
For convenience, rather than measuring the deviation from exact chiral symmetry using fh, 
which is scale dependent, we shall use m%- 

In terms of m n , current lattice calculations are typically restricted to pion masses larger 
than 500 MeV, with some pioneering work reporting preliminary results as low as 310 MeV. 
In order to compare these results with experimental data on hadron properties it is necessary 
to extrapolate the calculations at large pion masses to the physical value. In doing so it 
is crucial to respect the constraints imposed by chiral symmetry in QCD. In particular, as 
we discuss below, the existence of Goldstone bosons necessarily leads to behaviour which is 
non-analytic in the quark mass. 

The structure of this article is that we first explain the origin of the non-analyticity 
associated with Goldstone boson loops. We then explain, using the specific case of the 
nucleon mass, how this non-analytic structure has been incorporated into a new method for 
extrapolating hadron masses from the large values characteristic of lattice calculations to 
the physical region. The consequences of this for the sigma commutator are then explained. 
Next we turn to recent results for baryon electromagnetic properties. Finally we discuss 
the most recent investigations of the proton structure function, especially the importance 
of chiral symmetry in connecting existing calculations of lattice moments with data. We 
conclude with a summary of the promised insights into the nature of hadron structure within 
QCD that follow from all these investigations. 



For our purposes the primary significance of spontaneous chiral symmetry breaking in 
QCD is that there are contributions to hadron properties from loops involving the resulting 
Goldstone bosons. These loops have the unique property that they give rise to terms in an 
expansion of most hadronic properties as a function of quark mass which are not analytic. 
As a simple example we consider the nucleon mass. The most important chiral corrections 
to M N come from the processes N — > Ntt — > N (ojvjv) and iV — > Air — > N (ct/va)- (We 
will come to what it means to say these are the most important shortly.) We write = 
j^bare _|_ _|_ j n ^he heavy baryon limit one has 



where g^, f n are strictly evaluated in the chiral limit. Here u(k) is a natural high momentum 
cut-off which is the Fourier transform of the source of the pion field (e.g. in the cloudy bag 
model (CBM) it is 3ji (kR) / kR, with R the bag radius ) . From the point of view of PCAC 
it is natural to identify u(k) with the axial form factor of the nucleon, a dipole with mass 
parameter 1.02 ± 0.08GeV. 

Regardless of the form chosen for the ultra-violet cut-off, one finds that u^n is a non- 
analytic function of the quark mass. The leading non-analytic (LNA) piece of u^n is inde- 
pendent of the form factor and gives 



II. GOLDSTONE BOSON LOOPS AND NON-ANALYTICITY 
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This has a branch point, as a function of m, at m — 0. Such terms can only arise from 
Goldstone boson loops. 



It is natural to ask how significant this non-analytic behaviour is in practice. If the pion 
mass is given in GeV, <x^ A = — 5.6m^ and at the physical pion mass it is just -17 MeV. 
However, at only three times the physical pion mass, m n = 420MeV, it is -460MeV - half 
the mass of the nucleon. If one's aim is to extract physical nucleon properties from lattice 
QCD calculations this is extremely important. The most sophisticated lattice calculations 
with dynamical fermions are only just becoming feasible at such low masses and to connect 
to the physical world one must extrapolate from m n ~ 500MeV to m w = 140MeV. Clearly 
one must have control of the chiral behaviour. 

Figure [1] shows recent lattice calculations of Mn as a function of from CP-PACS and 
UKQCD Q]. The dashed line indicates a fit which naively respects the presence of a LNA 
term, 



with a, (3 and 7 fitted to the data. While this gives a very good fit to the data, the chiral 
coefficient 7 is only -0.761, compared with the value -5.60 required by chiral symmetry. If 
one insists that 7 be consistent with QCD the best fit one can obtain with this form is the 
dash-dot curve. This is clearly unacceptable. 

An alternative suggested recently by Leinweber et al. J5|, which also involves just three 
parameters, is to evaluate ctnn and ctjva with the same ultra-violet form factor, with mass 
parameter A, and to fit Mn as 



Using a sharp cut-off (u(k) = 9(A—k)) these authors were able to obtain analytic expressions 
for (Tjvtv and una which reveal the correct LNA behaviour - and next to leading (NLNA) in 
the Att case, ctna NA ~ Tn^hxm^. 



A. Case Study: the Nucleon Mass 



M N = a + Pml + 777^ 



(4) 



M N = a + f3ml + a NN {m % , A) + o-jva^tt, A). 



(5) 
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FIG. 1. A comparison between phenomenological fitting functions for the mass of the nucleon 
- from Ref. ||. The two parameter fit corresponds to using Eq.(||) with 7 set equal to the value 
known from xPT. The three parameter fit corresponds to letting 7 vary as an unconstrained fit 
parameter. The solid line is the two parameter fit based on the functional form of Eq.([|). 

These expressions also reveal a branch point at = Ma — Mjy, which is important if 
one is extrapolating from large values of m ff to the physical value. The solid curve in Fig. 1 
is a two parameter fit to the lattice data using Eq.(^), but fixing A at a value suggested by 
CBM simulations to be equivalent to the preferred 1 GeV dipole. A small increase in A is 
necessary to fit the lowest mass data point, at m 2 ~ 0.1 GeV 2 , but clearly one can describe 
the data very well while preserving the exact LNA and NLNA behaviour of QCD. 




B. Consequences for the Sigma Commutator 

The analysis of the lattice data for Mn, incorporating the correct non-analytic behaviour, 
can yield interesting new information concerning the sigma commutator of the nucleon: 

<tn = \{N\[Qx, [Qis,H QCD ]]\N) = (N\m(uu + dd)\N). (6) 

This is a direct measure of chiral SU(2) symmetry breaking in QCD, and the widely accepted 
experimental value is 45 ± 8MeV |J. (Although there are recent suggestions that it might 
be as much as 20 MeV larger 0.) Using the Feynman-Hellmann theorem one can also write 

dm 71 <9m 2 

Historically, lattice calculations have evaluated < N\{uu + dd)\N > at large quark mass and 
extrapolated this scale dependent quantity to the "physical" quark mass, which had to be 
determined in a separate calculation. The latest result with dynamical fermions, — 18 ±5 
MeV [§], illustrates how difficult this procedure is. On the other hand, if one has a fit to 
Mm as a function of which is consistent with chiral symmetry, one can evaluate <ttv 
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directly using Eq.(0). Using Eq.fl5|) with a sharp cut-off yields ~ 55 MeV, while a dipole 
form gives ~ 45 MeV ||. The residual model dependence can only be removed by 
more accurate lattice data at low rrv^. Nevertheless, the result G (45,55) MeV is in 
very good agreement with the data. In contrast, the simple cubic fit, with 7 inconsistent 
with chiral constraints, gives ~ 30 MeV. Until the experimental situation regarding ctjv 
improves, it is not possible to draw definite conclusions regarding the strangeness content of 
the nucleon. However, the fact that two-flavour QCD reproduces the current prefered value 
should certainly stimulate some thought and a lot of work. 



III. ELECTROMAGNETIC FORM FACTORS 
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FIG. 2. Recent data for the neutron electric form factor in comparison with CBM calculations 
for a confining radius around 0.95fm - from Ref. flCf| . 

It is a general consequence of quantum mechanics that the long-range charge structure 
of the proton comes from its 7r + cloud (p — ► mr + ), while for the neutron it comes from its 7T~ 
cloud (n — > p7i~). However, it is not often realized that the LNA contribution to the nucleon 
charge radius goes like lnm^ and diverges as rh — >• [IT]. This cannot be reproduced by a 
constituent quark model. Figure § shows the latest data from Mainz and NIKHEF for the 
neutron electric form factor, in comparison with CBM calculations for a confinement radius 
between 0.9 and 1.0 fm. The long-range 7r~ tail of the neutron plays a crucial role. 

While there are only limited (and indeed quite old) lattice data for hadron charge radii, 
recent experimental progress in the determination of hyperon charge radii has led us to 
examine the extrapolation procedure for obtaining charge data from the lattice simulations 
rr^] . Figure ^] shows the extrapolation of the lattice data for the charge radius of the 
proton. Clearly the agreement with experiment is much better once the chiral log required 
by chiral symmetry is correctly included, than if, for example, one simply made a linear 
extrapolation in the quark mass (or m%). Full details of the results for all the octet baryons 
may be found in Ref. [Q| . 
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FIG. 3. Fits to lattice results for the squared electric charge radius of the proton - from Ref. 
[12]. Fits to the contributions from individual quark flavours are also shown: the n-quark sector 
results are indicated by open triangles and the d-quark sector results by open squares. Physical 
values predicted by the fits are indicated at the physical pion mass, where the full circle denotes 
the result predicted from the first extrapolation procedure and the full square denotes the baryon 
radius reconstructed from the individual quark flavor extrapolations. (N.B. The latter values are 
actually so close as to be indistinguishable on the graph.) The experimental value is denoted by 
an asterisk. 



The situation for baryon magnetic moments is also very interesting. The LNA contribu- 
tion in this case arises from the diagram where the photon couples to the pion loop. As this 
involves two pion propagators the expansion of the proton and neutron moments is: 

^ = ^ n) T am n + 0(ml). (8) 

Here //{) is the value in the chiral limit and the linear term in m n is proportional to 
a branch point at fh — 0. The coefficient of the LNA term is a = 4.4/ijvGeV _1 . At the 
physical pion mass this LNA contribution is 0.6/ijv, which is almost a third of the neutron 
magnetic moment. No constituent quark model can or should get better agreement with data 
than this. 

Just as for M N} the chiral behaviour of (J?™ is vital to a correct extrapolation of lattice 
data. One can obtain a very satisfactory fit to some rather old data, which happens to be 
the best available, using the simple Pade ||14|| : 



p(n) 

uftri = ^2 (g) 
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FIG. 4. Absolute value of the ratio of the proton to neutron magnetic moments as a function 
of m 2 obtained from the Pade approximants in Eq. We stress that the behaviour as m\ — > 
is model independent. 

The data can only determine two parameters and Eq.(^) has just two free parameters 
while guaranteeing the correct LNA behaviour as — > and the correct behaviour of 
heavy quark effective theory (HQET) at large m 2 . The extrapolated values of \x v and [i n at 
the physical pion mass, 2.85 ± 0.22fi N and —1.90 ± 0.15/itv, respectively, are currently the 
best estimates from non-perturbative QCD [14]. For more details of this fit we refer to Ref. 



I4| , while for the application of similar ideas to other members of the nucleon octet we refer 
to Ref. jl5| , and for the strangeness magnetic moment of the nucleon we refer to Ref. jn| . 

Incidentally, from the point of view of the naive quark model it is interesting to plot 
the ratio of the absolute values of the proton and neutron magnetic moments as a function 
of m 2 . The agreement of the constituent quark result, namely 3/2, with the experimental 
value to within a few percent is usually taken as a major success. However, we see from Fig. 
H that it is in fact fortunate to obtain such close agreement JT7] . We stress that the large 
slope of the ratio near m 2 = is model independent. 



IV. STRUCTURE FUNCTIONS 



The parton distribution functions (PDFs) of the nucleon are light-cone correlation func- 
tions which, in the infinite momentum frame, are interpreted as probability distributions 
for finding specific partons (quarks, antiquarks, gluons) in the nucleon. They have been 
measured in a variety of high energy processes, ranging from deep-inelastic lepton scattering 
to Drell-Yan and massive vector boson production in hadron-hadron collisions. A wealth of 
experimental information now exists on spin-averaged PDFs, and an increasing amount of 



data is being accumulated on spin-dependent PDFs |L8 

At high momentum transfer (Q 2 ) the dominant component of the PDFs are determined 
by non-perturbative matrix elements of certain "leading twist" operators. In principle these 
matrix elements, which correspond to moments of the measured structure functions, con- 
tain vital information about the non-perturbative structure of the target. An extensive phe- 
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nomenology has been developed over the years within model QCD studies, and in some cases 
remarkable predictions have been made from the insight gained into the non-perturbative 
structure of the nucleon. An example is the d — u asymmetry, predicted [|llj on the basis of 
the nucleon 's pion cloud [p0| |, which has been spectacularly confirmed in recent experiments 
at CERN and Fermilab JH]]. Other predictions, such as asymmetries between strange and 
antistrange [22] and spin-dependent sea quark distributions, Au — Ad, still await experi- 
mental confirmation. Note that none of these could be anticipated without insight into the 
non-perturbative structure of QCD. 

Despite the phenomenological successes in correlating deep-inelastic and other high en- 
ergy data with low energy hadron structure, the ad hoc nature of some of the assumptions 
made in deriving the low energy models from QCD leaves open a number of questions about 
the ability to reliably assign systematic errors to the model predictions. One approach in 
which structure functions can be calculated systematically from first principles, and which 
at the same time allows one to search for and identify the relevant low energy QCD degrees 
of freedom, is lattice QCD. 
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FIG. 5. First moment of the difference u — d from various lattice QCD simulations (QCDSF 
[24-26] and MIT [27]), at a scale Q 2 = 4 GeV 2 . Calculations from the CBM are shown as small 
squares. The dashed curve is a simple fit which is linear in m 2 , while the solid curve incorporates 
the constraints of chiral symmetry, as in Eq.(lO). 



Early calculations of structure function moments within lattice QCD were performed 



by Martinelli and Sachrajda ||23|| . However, the most comprehensive analysis has been 
performed by the QCDSF Collaboration [24-26] - albeit within quenched QCD. Recently 
the MIT group has performed the first full (unquenched) QCD calculations of non-singlet 
moments |27| . The moments from the full QCD simulations are very similar to those from 
the quenched calculations. This is consistent with the suggestions of chiral quark models, 
like the CBM, that in the mass region currently accessible quark loops are suppressed. 

As for the other nucleon properties discussed above, we propose to extraplate the lattice 
data to the physical pion mass using a formula which is compatible with the LNA structure 



S 



of the PDFs. This behaviour was derived recently, with the result that the LNA behaviour 



involved a term in lnm x [28]. For an initial investigation we concentrate on the non- 
singlet combination of PDFs, u— d, in which "disconnected" quark loops cancel. Calculations 
based on the CBM (which incorporate the LNA chiral structure just discussed) actually 
produce quite a reasonable description of the behaviour of the moments of the PDFs as 
a function of quark mass, as shown in Fig. |5] (open squares). More important from the 
phenomenological point of view, the CBM calculations (for the n'th moment of the PDFs) 
can be fit with the simple expansion in m n : 



2 

m„ 



(xl - xf) =a n + b n m n + a n c hNA m n In - f , (10) 

where clna is model independent. 

The scale \i in Eq.(|T0"D is effectively the scale at which the rapid, chiral variation at 
low m n turns off. The best fit to the lattice data is obtained with a value fi ~ 0.4 — 0.5 
GeV - a very similar scale to that found, for example, for the magnetic moments. Clearly 
Eq . (|i~0"D gives a very good description of the lattice data for the first moment of the non- 
singlet distribution d — u. Taking into account the rapid chiral variation as m 2 — > there 
is also quite good agreement between the extrapolated value of the first moment and the 
experimentally determined moment. A similar result holds for the second and third moments 
too [El. 



V. CONCLUSION 

In the light of the numerous examples presented in this brief review, it should be evident 
that the study of hadron properties as a function of quark mass shows a clear pattern: 

• In the region of quark masses m > 60 MeV or so (m 7T greater than typically 400- 
500 MeV) hadron properties are smooth, slowly varying functions of something like a 
constituent quark mass, M ~ Mq + cm (with c ~ 1). 

• Indeed, Mn ~ 3M, M p ^ ~ 2M and magnetic moments behave like 1/M. 

• As m decreases below 60 MeV or so, chiral symmetry leads to rapid, non-analytic 
variation, with SMn ~ m 3//2 , S/jh ~ m 1 / 2 and 

5 < r 2 >ch~ him. 

• Chiral quark models like the cloudy bag provide a natural explanation of this transi- 
tion. The scale is basically set by the inverse size of the pion source - the inverse of 
the bag radius in the bag model. 

These are remarkable results that will have profound consequences for our further explo- 
ration of hadron structure within QCD as well as the analysis of the vast amount of data now 
being taken concerning unstable resonances. In terms of immediate results for the structure 
of the nucleon, we note that the careful incorporation of the correct chiral behaviour of QCD 
into the extrapolation of its properties calculated on the lattice has produced: 

• The best values of the proton and neutron magnetic moments from QCD. 
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• The best value of the sigma commutator. 



• Improved values for the charge radii of the baryon octet. 

• Improved values for the magnetic moments of the hyperons. 

• Good agreement between the extrapolated moments of the non-singlet distribution 
u — d and the experimentally measured moments. 

In addition, although we did not have time to discuss it, this approach has led to the best 
current value for the strangeness magnetic moment of the proton from lattice QCD 



Clearly, while much has been achieved, even more remains to be done. It is vital that 
lattice calculations with dynamical fermions are pushed to the lowest possible quark masses, 
taking advantage of developments of improved actions and so on. It is also vital to further 
develop our understanding of the physics of chiral extrapolation by comparison with these 
new calculations, by looking at new applications and by further comparison with chiral 
models. 
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